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ABSTRACT 

(N , 

The narrow emission line spectra of active galactic nuclei are not accurately 
^ _ described by simple photoionization models of single clouds. Recent Hubble 

' Space Telescope images of Seyfert 2 galaxies show that these objects are rich 

oo ■ 

O 



with ionization cones, knots, filaments, and strands of ionized gas. Here we 
lO ' extend to the narrow line region the "locally optimally emitting cloud" (LOG) 

^ ' model, in which the observed spectra are predominantly determined by powerful 

^ ■ selection effects. We present a large grid of photoionization models covering a 

■ wide range of physical conditions and show the optimal conditions for producing 

I ' many of the strongest emission lines. We show that the integrated narrow 

^ ■ line spectrum can be predicted by an integration of an ensemble of clouds, 

^ ■ and we present these results in the form of diagnostic line ratio diagrams 

making comparisons with observations. We also predict key diagnostic line 
i. ratios as a function of distance from the ionizing source, and compare these to 

j_j ■ observations. The predicted radial dependence of the [Qui] /[Oil] ratio may be 

^ matched to the observed one in NGG 4151, if the narrow line clouds see a more 

intense continuum than we see. The LOG scenario when coupled with a simple 
Keplerian gravitational velocity field will quite naturally predict the observed 
line width versus critical density relationship. The influence of dust within the 
ionized portion of the clouds is discussed and we show that the more neutral gas 
is likely to be dusty, although a high ionization dust-free region is most likely 
present too. This argues for a variety of NLR cloud origins. 

Subject headings: galaxies: Seyfert — line: format ion 
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1. Introduction 

The failure of "simple" photoionization models to describe the narrow emission line 
spectrum of active galactic nuclei (AGN) was shown dramatically by Filippenko & Halpern 
(1984), who found that the line widths of the optical forbidden lines were correlated 
with the critical density and ionization potential of those lines. This correlation indicates 
stratification of the photoionized region, with the lines emitting near their critical densities. 

Recent Hubble Space Telescope images of Seyfert 2 galaxies show that these objects are 
rich with ionization cones (Evans et al. 1991; Tsvetanov & Walsh 1992; Wilson et al. 1993; 
Macchetto et al. 1994; and Arribas et al. 1996), knots (Acosta-Pulido et al. 1996), and 
filaments, and strands (Falcke et al. 1996) of ionized gas. We should not expect that these 
complex environments can be described with simple, single-component photoionization 
models. Similar images of Seyfert 1 type objects show a less complex, symmetric structure 
that may indicate a viewing angle preference consistent with the unification model of 
Seyfert galaxies (Schmitt & Kinney 1996). 

Two-component photoionization models have been proposed by Binette, Wilson, & 
Storchi-Bergmann (1996) and have met with some success in describing the emission line 
spectrum of Seyfert galaxies. This model assumes that the narrow line region (NLR) is 
composed of both matter-bounded and ionization-bounded clouds, and that the covering 
fraction of the two populations can account for the range of observed line emission. A 
multi-component photoionization model has been proposed by Komossa & Schulz (1997), 
who assume that the NLR is effectively composed of a few clouds with a range of gas 
densities and distances from the central AGN. Both of these types of models do much 
better at reproducing the scatter in the line ratio observations than single component 
photoionization models. 

An alternative to photoionization by the central source is photoionization by fast 
moving shocks. Dopita & Sutherland (1995) describe in detail these "photoionizing" shocks 
and argue that all narrow emission lines can be excited by these processes. However, 
Morse, Raymond, & Wilson (1996) and Wilson (1997) discuss the viability of such shocks 
in describing the emission spectrum of active galaxies. In particular, Wilson (1997) argues 
that the energetics of radiative shocks may not account for the observed line luminosities. 
It is still unclear exactly what role shocks play in the NLR environment; that is, whether 
they provide the bulk of the energy input for the observed emission or merely compress gas 
so that it emits more strongly when photoionized by a central source (Pogge 1997). 

We will assume that photoionization by the central source is the dominating mechanism 
responsible for NLR emission. We extend the locally optimally-emitting clouds (LOG) 
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hypothesis of Baldwin et al. (1995) to NLR gas. This assumes that in nature, gas with 
a wide range of physical conditions is present and the emission line spectrum we observe 
is predominately a result of a simple, yet powerful, selection effect: we observe lines from 
clouds best able to emit them. To this end we present the results of a large grid of 
photoionization models and make predictions of the nature of the integrated spectrum. 

2. Model calculations 
2.1. Assumptions 

For simplicity we will assume that the source of narrow line emission observed in active 
galaxies is from photoionized clouds that are similar to illuminated molecular clouds in our 
own galaxy. The source of these clouds may be four-fold: (1) molecular cloud complexes 
in the disk of the galaxy; (2) the dusty molecular torus that is beheved to surround the 
central source and which will certainly be photoionized; (3) ISM entrained behind a radio 
jet so that photoionization by the central source becomes important; and (4) material 
"snow plowed" in the form of bow shocks in front of a radio jet such that it too may be 
photoionized by the central source. We will ignore emission that may be due to shock 
excitation coming from photo-ionizing shocks, although how important this emission is to 
the overall spectrum is unclear (Wilson 1997). 

We have used the spectral synthesis code CLOUDY (version 90.03; Ferland 1996) to 
calculate the emission from plane parallel, constant hydrogen density clouds ionized by a 
continuum similar to that expected from a typical Scyfcrt galaxy with Lion = lO'*'^ ''"^ ergs s~^. 
The shape of the ionizing continuum was chosen to be a combination of a UV-bump of the 
form fi, oc v~^-^exp{—hu / kTcut) and an X-ray power law of the form f^, oc ° spanning 
13.6 eV to 100 keV. The UV-bump cutoff temperature, Tc„t, was chosen such that the 
UV-bump peaked (in vF^) at 48 eV. The UV and X-ray components were combined with a 
typical Seyfert UV to X-ray spectral slope, ctox = —1.2. 

The ionization/thermal equilibrium and radiative transfer calculations for a single 
cloud proceeded until one of the following three conditions were met. First, the electron 
temperature dropped below 3000 K. Gas with lower temperatures does not contribute 
significantly to the emission lines presented in this paper. Second, we do not allow clouds 
to have thicknesses in excess of 10% of its distance from the central continuum source. 
Third, we do not allow the total hydrogen column density to exceed 10^"^ cm~^; this mainly 
affects the very high ionization clouds and prevents them from becoming Thomson thick. 
In practice the third condition had little impact on the narrow emission lines. The largest 
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cloud of the grid presented below is one with a thickness of ~16 pc, consistent with the size 
of large molecular clouds. We will comment further on the effects of these stopping criteria 
in a later section. 

Finally, we assume that each cloud sees the full continuum with no obscuration. 
At large enough distances, intervening clouds or diffuse ISM may attenuate the ionizing 
spectrum significantly. We have looked at the effects of such ISM attenuation on the 
incident continuum: An ISM with gas densities as low as 0.1 cm~^ such as found in the local 
ISM of the Milky Way (Wood & Linsky 1997) would have no effect on the results presented 
below. A higher density ISM (n(II) = 0.5 cm~^) however, would have an effect at large 
distances from the ionizing source, due to the larger column density and lower ionization 
resulting in a higher gas opacity. For the higher density gas and the shape and luminosity 
of the incident continuum given above, the attenuated continuum has an optical depth of 
unity at 4 Ryd at a distance of logR ~ 21.75 cm from the ionizing source, resulting in 5 - 
20% changes in the line luminosities. At larger distances the attenuation by a dense ISM 
would result in a decrease in the luminosity of the emitted spectrum of up to a factor of 10. 
We neglect this complication for simplicity. 



2.2. Reprocessing efRciency 

We have computed, on the University of Kentucky Convex Exemplar parallel 
supercomputer, a grid of 1,881 of the single cloud models described above in a plane of 
hydrogen gas density and distance from the central ionizing source. We chose a range of 
parameters in order to cover all the possible physical conditions of narrow line region gas. 
The hydrogen gas density, n{H), ranged from 10^ - 10^° cm^^. The incident continuum 
flux was varied so that for an assumed ionizing continuum luminosity Lion = lO''^'^ ergs s~^ 
(monochromatic luminosity logL;^(1450 A) = 39.558 ergs s~^ for the continuum shape 
used here) the distance from the ionizing source of radiation, R, varied from 10^^ - 10^^ cm. 

1/2 

This distance range scales with the luminosity as . The conversion from incident flux to 
luminosity and radius assumes that the ionizing radiation is emitted isotropically. However, 

observations of ionizing "cones" (Evans et al. 1991; Tsvetanov & Walsh 1992; Wilson et 
al. 1993; Macchctto et al. 1994; and Arribas et al. 1996) and UV-photon deficits (Binette, 
Fosbury, & Parker 1993; and Morse, Raymond & Wilson 1996) imply that the gas may see 
a different ionizing continuum than we do. 

Calculations were performed assuming two abundance sets: solar and a set including 
dust that varied with distance from the ionizing source. We will describe each set of 
abundances in turn and compare the two in terms of diagnostic line ratio diagrams in a 
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later section. 



2.2.1. Solar abundances 

Figure 1 sliows tiie results of grid calculations, for solar abundance clouds, in the form 
of contour plots of logarithm equivalent width referred to the incident continuum at 4860 A. 
The solar abundances are from Grevesse & Anders (1989) and Grevesse & Noels (1993) and 
have the following values: 

H :1.00E+00 He:1.00E-01 Li:2.04E-09 Be:2.63E-ll B :7.59E-10 C :3.55E-04 N :9.33E-05 
O :7.41E-04 F :3.02E-08 Ne:1.17E-04 Na:2.06E-06 Mg:3.80E-05 Al:2.95E-06 Si:3.55E-05 
P :3.73E-07 S :1.62E-05 Cl:1.88E-07 Ar:3.98E-06 K :1.35E-07 Ca:2.29E-06 Sc:1.58E-09 
Ti:1.10E-07 V :1.05E-08 Cr:4.84E-07 Mn:3.42E-07 Fe:3.24E-05 Co:8.32E-08 Ni:1.76E-06 
Cu:1.87E-08 Zn:4.52E-08 

Furthermore, we assume that there are no dust present in the solar abundance computations. 

The ionization parameter provides a useful homology relationship with the emission 
lines (Davidson 1977). U{H) is defined as the ratio of hydrogen ionizing photon density 
to hydrogen density, U{H) = ^{H)/n{H)c, where ^{H) is the flux of hydrogen ionizing 
photons and c is the speed of hght. For reference a dashed hne is placed on the upper left 
hand plot of the first panel representing \ogU{H) — 3.0. The ionization parameter increases 
from top right to bottom left in these diagrams; lines of constant U{H) have a slope of 
dlog(R)/dlog(n) = -0.5. 

Shown in Figure 1 are 23 strong optical, UV and infrared recombination, resonance, 
and forbidden hues. The temperature of the illuminated face of the cloud is also shown 
in the upper right hand plot of the first panel. For many emission lines a ridge of near 
maximum equivalent width runs diagonally across the distance - gas density plane, roughly 
parallel to lines of constant ionization parameter. Most lines are emitted optimally for a 
narrow range of ionization parameter spanning this ridge. For larger (smaller) ionization 
parameters the gas is over (under) ionized, and the line is not efficiently emitted. This 
explains the sudden drops in the emission line equivalent widths on either side of their 
ridges. Moving along the ridge to increasing gas densities, at near constant ionization 
parameter, the forbidden lines become coUisionally deexcited and the line equivalent width 
falls off. Moving along the ridge at constant ionization parameter to smaller gas densities, 
the equivalent widths of some lines also diminish as other lines of similar ionization, but 
lower critical density, become important coolants. This effect is especially dramatic for C IV 
A1549 and Ovi A1035. 
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The stopping criteria assumed in §2.1 do affect the resulting spectrum. The grid of 
models shown in Figure 1 includes clouds that are both optically thin and optically thick to 
hydrogen ionizing radiation. Optically thin clouds generally occupy the lower left portion 
of the distance - density plane below a line with slope dlog(R)/dlog(n) ~ —2/3 and 
intersecting the lower right hand corner. The temperature of this gas (see Fig. 1) ranges 
from 20,000 K to several million K. Emission lines presented in this paper are not emitted 
by gas with ^ 10^ K (X-ray lines are). These hot clouds are generally truncated by 
the plane parallel condition discussed earlier. Optically thick clouds are in the upper right 
portion of the figures, and the gas generally has nebular temperatures, except the extreme 
upper right corner where the electron temperatures can be as low as a few hundred Kelvin 
and the gas is mostly molecular. Generally the optically thick clouds were stopped because 
the electron temperature fell below 3000 K (§2.1). The total hydrogen column density 
restriction mentioned in §2.1 applies to a small number of high density optically thick clouds 
in a triangular region bounded by a distance of 1 pc on the the topside and a line roughly 
parallel to log U ^ — 1.0 on the leftside. The clouds in this region do not significantly 
contribute to the narrow forbidden line spectrum. 

Different types of lines can be seen in Figure 1 to be emitted from clouds with a range 
in gas densities and distances. Recombination lines, such as HP and La, generally form 
broad planes of maximum equivalent width whereas the forbidden lines, such as [Ne v] , are 
either narrow ridges or peaked islands ([On] A7325). Figure lb shows the differences in the 
nebular and auroral lines of [O ll] and [O ill] . Notice that the auroral features are formed 
optimally at higher densities, because of their higher critical density. Figure Ic shows the 
contour plots of a few of the stronger UV-lines. Figure Id shows a few strong [Fell] infrared 
fine-structure lines from a 16-lcvel model atom. The model atom has the same configuration 
as that shown in Thompson (1995) and uses the collision data of Pradhan & Zhang (1993) 
and Einstein As from Quinet, Lc Dourneuf & Zeippen (1996). Vcrncr ct al. (1997) compare 
the differences between the 16- level model atom and a much larger definitive 371-level 
calculation. Preliminary comparisons show that predicted line intensities of the 16-level 
atom differ typically by a few percent and at maximum by 10% from the larger 371-level 
calculation. 

To illustrate the problems encountered when trying to use conventional diagnostic 
line intensity ratios for an ensemble of clouds with a wide range of physical conditions, we 
momentarily assume that each emission line in the total spectrum is produced with the 
maximum possible equivalent width found in the figures. We then use this spectrum to 
reproduce the standard density and temperature diagnostic line ratios. This is the simplest 
way to document the effects of a distribution of clouds. Three density indicators are 
given by Osterbrock (1989). These include the [Sll] A6716/A6731 and [Oil] A3729/A3726 
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doublets and the [Cm] A1907/ Cm] A1909 ratio. For [Fell], Pradhan & Zhang (1993) use 
the 1.534 /im/1.644 //m ratio as a density indicator. Assuming a temperature of 10,000 K, 
the [On] ratio we predict gives an Ue of 300 cm~^, for the predicted [Sll] doublet the 
density is 700 cm~^, for the [Fell] ratio Ug is 20,000 cm^^, and the [Cm]/ Cm] ratio gives a 
density of nearly 300,000 cm'^ The [O III] (A4959 + A5007)/A4363 temperature - density 
diagnostic using the peak emission of Figure 1 and an electron temperature of 10,000 K 
gives an electron density on the order of a few million cm~^. The point of this exercise is 
clear: line diagnostics fail when the gas has a wide range of physical conditions. As Figure 1 
shows, many of these lines form at different densities and temperatures. However, by using 
the diagnostic line ratios, one assumes that the lines are formed in gas of similar physical 
conditions, which is obviously not the case in Figure 1, and is likely not the case in nature. 
We will in a later section integrate the "visibility functions" of Figure 1 to calculate an 
emission spectrum by assuming that the observed spectrum is the result of an ensemble of 
clouds with differing properties. 

2.2.2. A dusty NLR 

In the previous section we have assumed that the narrow line region is a dust free 
environment. We now recompute the grid assuming that the clouds arc dusty. Grains will 
not survive in an environment with high photoionizing flux due to the various sublimation 
processes discussed by Laor & Drainc (1993). Therefore, after the manner of Netzer & 
Laor (1993) we vary the abundances as a function of distance from the ionizing source, 
in order to mimic these destruction processes. At small distances from the central source 
we assume the abundances are solar, as given in the previous section. At a distance of 
approximately 10^^'^ cm from the central source (for the ionizing luminosity given above), 
Orion type graphite grains (Baldwin et al. 1991) are just at their sublimation temperature 
(~1750 K) at the cloud face. At distances larger than this we include graphite grains within 
the calculation of the single cloud and deplete carbon such that its abundance is typical of 
an Hll region (see below and Baldwin et al. 1996). Orion type silicate grains sublimate at 
~1400 K, corresponding to a distance of lO^'''^ cm. Beyond this radius the abundances are 
set to values which approximate those found in the Orion nebula: 

H :1.00E+00 He:9.50E-02 C :3.00E-04 N :7.00E-05 O :4.00E-04 Ne:6.00E-05 Na:3.00E-07 
Mg:3.00E-06 Al:2.00E-07 Si:4.00E-06 S :1.00E-05 Cl:1.00E-07 Ar:3.00E-06 Ca:2.00E-08 
Fe:3.00E-06 Ni:1.00E-07 

These abundances are based upon the results of several recent studies of the Orion nebula 
(Baldwin et al. 1991, Rubin et al. 1991, 1992, and Osterbrock, Tran & Veilleux 1992). 
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Although this discontinuous "turning on" of grains is unphysical and our choices of grain 
composition and gas depletions are specific, it serves as an example to simulate the presence 
of dust in the narrow line region. The grain physics used here is described by Baldwin et al. 



Figure 2 shows the same contour plots as Figure 1 except that the individual models 
include the grains as specified above. The discontinuous grain sublimation distances are 
easily seen in the figure as sharp discontinuities in the contour plots at logR — 16.9 and 



Comparisons of Figures 1 and 2 show four major effects of grains: (1) the depletion of 
refractory elements; (2) the weakening of emission lines due to absorption of the incident 
continuum by dust at large U{H); (3) the photoelectric heating of the gas by grains; and (4) 
in the case of resonance lines such as La A1216 and C IV A1549, line destruction by grains. 

Higher ionization lines, such as [Olll], [Nelll], and [Nev], tend to be emitted less 
efficiently than in the dust-free models, with the equivalent width falling by as much as 
a factor of 2. The weakening of these lines is caused by the absorption of the incident 
continuum by the dust. The lower ionization lines ([Sll], [Oil], [Nil], [O l] etc.) actually 
have their peak emission increase, though at lower ionization parameters, due to the 
photoelectric heating of the gas by the grains. Resonance line destruction is readily 
apparent in Figure 2c. The heavy depletion of iron is also apparent in Figure 2d, wherein 
the peak equivalent width drops by nearly a factor of ten. This last point will be discussed 
more fully in a later section. 



The equivalent widths shown in Figures 1 and 2 reflect the efficiency with which the 
clouds reprocess the incident continuum. It is clear that different lines are formed optimally 
in different places in the gas density - distance plane. If the distribution function of clouds 
in this plane is known, then the total line luminosity emitted by a set of clouds will be 
given by 



where F(r, n) is the emission line flux of a single cloud at radius r, gas density n, and 
ip{r,n) is the cloud distribution function, which is not necessarily analytical. 

The remaining question is the spatial distribution of the NLR clouds. Many type 
2 Seyferts have resolved NLR's, often with complex structures that are associated with 
linear or double-lobed radio sources (Bower et al. 1995; Pogge & De Robertis 1995; Capetti 



(1991). 



17.6 cm. 



3. Integrated line emission 




(1) 
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et al. 1996; Cooke et al. 1997), whereas the structure of the NLR in type 1 Seyferts is 
generally more compact and axisymmetric, at least in projection (Schmitt & Kinney 1996; 
Nelson et al. 1996). Since it is not clear how the emitting gas is distributed, in all its 
dimensions, we will assume that the distribution function il){r,n) can be approximated by 
ilj{r,n) oc f{r)g{n), where /(r) and g{n) are the cloud covering fractions with distance 
and gas density, respectively (see Baldwin et al. 1995). For simplicity we assume simple 
power laws: 

/(r) oc r'^ , g{n) oc n^. (2) 

The power law indices, 7 and /?, will be the only free parameters of our predicted integrated 
spectra. In the case where both weighting functions are proportional to the inverse of n or 
r, then the line equivalent widths of the clouds are equally weighted (in log space) across 
the gas density - distance plane. 

In the integrations computed below we assumed certain integration limits in the gas 
density - distance plane in order to exclude gas that may not be present in the NLR. 
We exclude gas with densities log n{H) > 8, the rough maximum critical density of the 
optical forbidden lines observed. By including gas with \ogn{H) ~ 9, the hue ratios do 
not change significantly except for those integrations that are weighted to higher densities 
(flatter g{n)). If present, the contribution of high density gas in the inner regions of the 
NLR cannot be significant or else the narrow line spectrum would appear vastly different 
(forbidden lines would be too weak relative to the hydrogen lines). We exclude gas that is 
closer than 0.1 pc to the ionizing source. This distance corresponds to the grain sublimation 
point described in §2.3. Integration limits were also included in ionization parameter space 
in order to exclude gas that is far over (under) ionized and does not contribute to the 
observed emission. These limits effectively exclude gas from the broad emission line region 
and very high density gas hundreds of parsecs away from the central engine. 

3.1. Diagnostic line ratio diagrams 

Integrations described above have been computed and plotted after the manner of 
Baldwin, Phillips & Terlevich (1981; hereafter BPT). Intensity ratios for the integrated 
model spectra are plotted in this manner in Figure 3. Two types of lines are shown in the 
figure. Solid fines represent dusty models with depleted refractory elements typical of H ll 
region type abundances (hereafter "dusty"). The dashed lines represent dust-free clouds 
with solar mctallicities (hereafter "dust-free"). Different values of g{n) (cf. Eq. [2]) are 
represented by different lines, with the power of f(r) ranging from —2.0 < 7 < 1.0 in 0.25 
increments along each line with the direction of more negative /(r) shown by the arrows. 



-10- 



For the dusty simulations (solid lines) three values of the g{n) index are shown (/3 = —1.0, 

— 1.6, —1.8), and the trend of more negative f3 are indicated by the appropriate arrows in 
the figures. For the dust-free (dashed lines) simulations the g{n) index are f3 = —1.0, —1.4, 

— 1.8, similarly the trends for the more negative index are indicated. The "S" and "L" 
indicate in the figures the average values of the line ratios for Seyfert, taken from Ferland & 
Osterbrock (1986), and LINER type spectra taken from Netzer (1990), with the exceptions 
of HellA4686, [Olll] A4363, from Ho, Filippenko & Sargent (1995) for the LINER galaxy 
M81. The observed range in [Oil] A7325 and [Slll] AA9069,9532 strengths were taken from 
Dopita & Sutherland (1995). 

To make it possible to follow particular models from panel to panel, we have marked 
two fiducial models in Figure 3. The stars represent a dusty model (sohd hues) with 
/(r) oc r~^'^ and g{n) oc n~^'^, while the squares represent a dust-free model (dashed lines) 
with /(r) oc r~^"^^ and g{n) oc n~^'^. These particular models were chosen because they 
approximately reproduce the [Oil] A3727/[Olll] A5007 and [Olll] A5007/H/3 hue ratios in 
the average Seyfert 2 spectrum. 

Interestingly, solar abundance, grainless calculations are not excluded by the line ratios 
presented in Figure 3. The dashed, dust-free, lines in Figure 3 tend to converge in the area 
of Seyfert observations. The dusty integrations tend to be shifted toward and into the 
LINER region of the diagrams, that is with higher [N ll] , [S ll] , and [O l] relative to Ha and 
higher [O ll] relative to [O ill] . 

Table 1 compares the fiducial models represented by the star and square in Figure 3 
with the observed average Seyfert 2 spectra from Ferland & Osterbrock (1986), Dopita & 
Sutherland (1995), Thompson (1995), and Villar-Martfn & Bincttc (1997). The UV/optical 
line ratios given by Ferland & Osterbrock, likely suffer from inaccurate absolute optical 
spectral flux calibrations coupled with substantially different aperture sizes between the 
lUE UV spectra and ground based optical spectra coupled with the extendcdness of the 
NLR, making many of these UV/optical narrow line ratios highly uncertain. The table also 
lists the observed range of equivalent width of H/3 taken from Netzer (1990) and Binette, 
Fosbury, & Parker (1993). 

In column 4 of Table 1 we give the dereddened UV/optical HST/FOS spectrum of the 
Seyfert 2 galaxy NGC 3393, from Cooke et al. (1997). The Cooke et al. observations were 
taken at two overlapping but significantly different aperture positions within a complicated 
two-dimensional hne-emitting structure; spectra covering 1100 A - 2300 A were taken at 
one position while the range 2300 A - 6800 A was covered at the other position. Because 
of the unintended positional offset it is impossible to unambiguously join together the two 
spectral regions. We have normalized the two sections of the NGC 3393 spectrum so that 
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the UV/optical Hell ratios would agree with case B predictions. This does not guarantee 
consistency amongst the other emission hues, but it is not clear that the uncertainty due 
in this case to a 0.4 arcsec offset between two 1 arcsec diameter apertures is any worse 
than the errors incurred through the comparison of lUE to ground based data, whereas the 
signahnoise in the NGC 3393 spectra is much higher than in the earlier data. 

Observations by Oliva et al. (1994) of the Seyfert type 2 Circinus galaxy are given in 
column 5 of Tabic 1. These observations include the ground-based optical and infrared 
spectra of many high ionization "coronal lines" (ionization potentials greater than 100 eV). 

Columns 6 and 7 of Table 1 show the results of the predicted integrated line spectrum 
for the dust-free and the dusty fiducial points. Recall that the integration parameters 
for the two fiducial models were chosen to simultaneously reproduce the [O III] /[Oil] and 
[OlIl]/H/3 emission line ratios of the average Seyfert galaxy, and with this in mind, the 
integrated spectrum is representative of the observed values. 

3.1.1. UV lines 

Compared to the mean Seyfert 2 spectrum the UV lines. La, CiV, and Cm], are 
generally undcrprcdictcd by the solar, dust-free integration shown in column 5 of Table 1. 
In column 6 of the table, the dusty integration underpredicts La, and CiV even worse than 
the dust-free case (generally because of resonance line destruction by dust grains). The 
exceptions being that Cm] A1909 fine, is within the observed range, and Mgll is on the low 
end of the range. The predictions of the weaker UV emission fines observed in NGC 3393, 
not included in the mean Seyfert 2 spectrum of Table 1, are generally consistent with the 
observations. The exceptions N v and He ll will be elaborated upon in the next section. 

The integrations of the dust- free clouds predict La/H/? ~28, which is close to the 
high density case B limit (~34), while the observations compiled by Ferland & Osterbrock 
(1986) suggest much larger ratios (30 - 70). Our calculations find a maximum Lq;/H/3 ~100 
occuring at \ogn{H) — 9.75 cm~^. Lq;/H/3 significantly larger than ~34 impfies the 

existance of a substantial population of high density clouds; more specifically 40 < La/H/3 
< 70 occurs for 5.5 < \ogn{H) < 8.5. However, the spectrum resulting from integrations 
that heavily weight this dense gas would not be consistent with the strong optical forbidden 
line strengths relative to H/9. 

Our predicted hydrogen line spectrum is in better agreement with the HST data of 
Cooke et al. (1997, see Table I). The corrected observed La/H/5 ~19 is inconsistent with 
the compilations of Ferland & Osterbrock (1986), with the caveat that this observed ratio 
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is somewhat uncertain. Lq;/H/3 ratios smaller than 23 are not possible, unless La is weak 
either due to very dense (> 10^^ cm~^) or dusty gas. The NGC 3393 Mgll/H/9 ratio is 
also consistent with the dusty simulations. HST observations of the narrow lines of several 
high luminosity AGN (Wills et al. 1993) also yield ratios that are more consistent with the 
dusty simulations. Clearly, uniform, long slit, high signal-to-noise, UV to infrared spectra 
are required to understand the narrow emission line spectra of AGNs. Future observations 
using the STIS on HST will accomplish a good part of this for the first time. 

3.1.2. Optical lines 

The predicted narrow optical forbidden lines are generally within the observed ranges 
of the mean spectrum with the exceptions of [Nelll] and [Nev] both being slightly 
underpredicted by ~50% and ~20%, respectively, indicating possible differences in 
abundances. However, these lines are predicted to be stronger relative to H/3 than observed 
in NGC 3393. The [S ll] , [S ill] lines are within the observed ranges of the mean Seyfert 2 
spectrum. 

Both integrations of the predicted nitrogen line [Nl], are consistent with the observed 
mean Seyfert 2 range, although the [Nil] line is predicted to be lower that observed by 
nearly half. Comparisons between both integrations and NGC 3393 of the nitrogen emission 
lines ([Nl], [Nil], and Nv) relative to H/3 indicate that nitrogen may be enhanced in that 
object relative to solar by a factor of ~2. 

The strength of the He ll A4686 permitted hne relative to H/3 is often used to constrain 
the ionizing continuum. A stronger Hell A4686 line (relative to H/3) is indicative of more 
He"*" ionizing photons. Simple photoionization models often underpredict the line while 
doing well with the rest of the optical spectrum (Ferland & Osterbrock 1986). Figure 3 and 
Table 1 show that the integrated model can predict this line within the observed range. 
Our fiducial integrations given in columns 5 and 6 of the table overpredict the entire He II 
spectrum relative to H/5 when compared to NGC 3393, indicative that NGC 3393 has a 
slightly softer continuum than the one assumed here (§2.1). 

The predictions of the relative strengths of the optical iron coronal lines, [FeVIl], [Fex], 
[FeXl], and [Fexiv] to H/? are in accord with the observations of the Circinus galaxy (Oliva 
et al. 1994) for the dust-free integrations. The dusty integrations underpredict these lines 
by factors of ~ 5 - 10. However, Oliva et al. (1994), Morisset & Pequignot (1996), and 
Ferguson et al. (1997) discuss the apparent large uncertainties in the coUision strengths 
of [Fex] A6375, [Fexi] A7892, and [Fexiv] A5303. Given these uncertainties we set the 
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coUision strengths for these transitions equal to 1 prior to running these simulations (note: 
Ferguson ct al. used the collision strengths from the literature). These lines will not be 
useful probes of the high ionization narrow line region until this issue is resolved. 

The weakness of the optical doublet [Call] AA7291,7324 has been cited as evidence 
for the presence of dust grains with the NLR (Ferland 1993; and Kingdon, Ferland & 
Feibelman, 1995). Tabic 1 indicates that the computations assuming solar abundances 
overpredict the [Call] A7291 line by factors of 3 - 160, with the dusty calculations on the 
very low end of the observations. This suggests that calcium is actually depleted by factors 
up to 160 compared to its solar abundance. 

3.1.3. Infrared lines 

The infrared forbidden [Fc ll] lines from the integrations described above are compared 
with the observations of NGC 4151 made by Thompson (1995). Because Thompson's 
observations included P/3, we can reference the [Fell] to H/9, as in Table 1, by assuming 
case B conditions for hydrogen. This is a valid assumption since P/3/H/9 predicted by our 
integrations is very nearly the case B ratio. Table 1 shows that the predicted intensities 
of the "dust-free" integration (column 3) are clearly too strong to match the observations, 
with some hues overpredicted by a factor of 10. Even for the "dusty" integration (column 
4), where the iron abundance is depleted by a factor of ten, relative to solar, the strongest 
[Fell] lines are too strong by a factor of 2 compared with H/3. This evidence suggests that 
iron may be depleted by factors up to 20 in NGC 4151. 

Simpson et al. (1996) and Veilleux et al. (1997) present IR spectroscopy of a number 
of Seyfert 2 galaxies. Simpson et al. found [Fell] 1.257 //m/P/3 1.21 (cr 0.81), and 

Veilleux et al. found [Fell] 1.257 /im/P/5 ^ 1.24 {a ^ 0.90). As in NGC 4151 these 
observed values are better matched to the predicted ratio from the dusty integration in 
Table 1 (~ 1.4). The dust-free integration predicts a ratio of ~ 7, much like the ratios 
arising in the shocked (grain free) environments of supernova remnants, plotted by both 
Simpson et al. and Veilleux et al. 

Although the [Fell] line ratios are too strong compared with the hydrogen lines, 
indicating an abundance difference, the relative [Fell] spectrum is fairly consistent with 
the observations. Thompson (1995) discusses at length the probable difficulties with the 
observation of weaker lines, such as 1.279 fim and 1.677 /im. 

Table 1 also includes predictions of the relative strength of the far-lR [Fell] 25.98/im 
line, the lowest transition within the ground term. This line is predicted to be quite 
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strong in the solar integration (4 times P/3). Recall that the simulations presented here 
excluded gas with Tg < 3000 K, and emission from these low excitation lines might still 
be important at lower electron temperatures. In fact, we find that this line and others 
that arise within the ground term will have substantial contributions from AGN-heated 
photodissociation regions, if present. None of the ~ 1 /xm [Fell] lines are emissive in gas 
temperatures below 3000 K. 

Our simulations of the infrared coronal lines shown in Table 1 are generally in accord 
with the observations of the Circinus galaxy for the dust-free integration. The [Sivi], 
[SiVIl], and [Siix] spectrum is well matched by the fiducial dust-free integrations. The 
dusty models do a poor job at matching the observations, particularly the [CaVlll] line. 
Our dust-free fiducial model predicts this line to within a factor of two, whereas the dusty 
integration underpredicts this fine by over a factor of 30. The predicted strength of the 
optical coronal fine [Fevil] A6087 is also consistent with the proposition that the high 
ionization lines are emitted in a dust-free environment. 



4. Line width and critical density 

Filippenko & Halpern (1984) showed that the fine width of optical forbidden lines 
is correlated with the critical density of those fines. Other observers (Fifippenko 1985, 
Appenzeller & Ostereicher 1988, Espey et al. 1994 and Ho, Filippenko, & Sargent 1996) 
have extended this work. It was realized by Filippenko & Halpern that the probable 
explanation for the observed correlation is that the lines emit most efficiently near their 
critical densities and that the line emission is spatially stratified. Our LOG models are a 
quantitative formulation of this idea. By just glancing at Figure 1, one observes that the 
peak emission of the set of lines shown occurs over a range of gas density and distance. If 
the dominant velocity field is one which diminishes with distance from the nucleus, then 
we expect that those lines that are emitted optimally at smaller distances (larger n{H)) 
to have broader line widths than those that are optimally emitted farther away (smaller 
n{H)) and qualitatively agree with the observed correlation. Consider the contour plots 
of [Olll] A4363 and [Sll] A6720 shown in Figure 1. It is clear by looking at their places of 
peak emission that [Olll] will be much broader than [Sll], as is observed. 

Figure 4 shows the results of a simple calculation of relative line widths assuming that 
the density - distance phase space system of clouds obey the covering fraction and density 
indices of the fiducial square shown in Figure 3. In lieu of a more realistic velocity field 
which takes into account the mass distribution of the bulge in addition to the supermassive 
black hole, we assume virial velocities about a simple point mass of 10^ Mq. In the figure 
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we show the full width at half maximum (FWHM) and the full width at 10% intensity 
(FW10%) for the ions indicated versus the logarithm of the critical density {nf^^). The 
slopes of the fitted linear regression in the figure were computed to be 0.193 (FW10%) and 
0.118 (FWHM), similar to those observed by Filippenko & Halpern (1984) and Espey at al. 
(1994), which range from 0.095 to 0.200. 

5. Application to Resolved NLR's. 

In a number of nearby Seyfert galaxies the NLR is spatially resolved even on 
ground-based long-slit and Imaging Fabry Perot spectra. These NLRs frequently have the 
form of either ionization cones (cf. Pogge 1988; Tadhunter & Tsvetanov 1989; Wilson 
et al. 1993; Robinson et al. 1994) or of apparent bow-shocks and cocoon shocks around 
linear radio sources (Bower et al. 1995; Pogge & De Robertis 1995; Capetti et al. 1996; 
Cooke et al. 1997). Typical scales are 200-400 pc arcsec"^ where h =Ho/100 km s^^ 
Mpc~^. With the imminent commisioning of the HST STIS spectrograph and ground-based 
adaptive optics systems, it should be possible to study spectroscopically a reasonable 
sample of NLR's with ~ 20 — 40 pc (~ 10^° cm) resolution. 

The exact scahng of such dimensions onto the y-axis of Figures 1 and 2 is comphcated 
by the fact that the (probably beamed) ionizing continuum luminosity seen by these NLRs 
is not directly measured. However, simple arguments based on counting recombination 
photons (cf. Binette, Fosbury & Parker 1993) imply that the gas sees an ionizing luminosity 
that is on average a factor of 10 larger than what we see. Therefore, the range in incident 
flux used in the model grid presented here is applicable to these resolved NLRs, and it is 
of interest to consider the predicted radial dependence of various line intensity ratios when 
integrated over the cloud gas density distribution. 

Figure 5 shows four different line ratios indicating the gas density, temperature, and 
ionization as functions of distance and gas density in the form of logarithmic contour plots. 
The figure also plots the ratios of the line intensities integrated over clouds of all densities, 
as functions of distance from the ionizing source; shown are the results of two different 
integrations along the density axis. 

5.1. Density 

The density indicator [Sll] A6716/A6731 is shown in the top two panels of Figure 5a. 
The line ratio has the linear value of 1.36 in the extreme upper left corner of the plot at 
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{logn{H) — 2.0, log = 22.0). Moving down and to the right along a hne of constant 
ionization parameter the solid lines are 1 dex increments and the dashed lines are 0.1 dex 
steps. The large flat plain of constant ratio corresponds to the high density limit or a linear 
value of ~ 0.44. 

The line ratio as a function of distance shows that for the integrations presented the 
high density limit is predicted at small distances. The solid line (integrated assuming that 
g{n) oc n"^") remains at that limit for larger distances, whereas the dashed line (integrated 
assuming a steeper density distribution, i.e. g{n) oc n^^'^) turns upwards sooner. The 
dashed curve reaches a peak vahic nearly equal to the low density limit given by Osterbrock 
(1989). The gas density inferred from the dashed curve falls off as R^^^^, while the density 
inferred from the solid curve is pinned at the high-density limit for larger values of R, then 
falls off more steeply than R~'^^^. 

5.2. Temperature 

The temperature indicator [Olll] A5007/A4363 (note that A4959 is not included) 
is shown in the bottom two plots of Figure 5a. The contour plot has a peak value of 
112 in the upper left hand corner (logn(^r) = 2.0, log i? = 21.0). Following a line of 
constant ionization parameter down and to the right, the solid contour line (which has a 
backwards "S" shape) has the logarithmic value of 2.0, and the dashed lines are 0.2 dex 
increments. The outer contour in the lower right part of the plot has a logarithmic value 
of —1.0 at {logn{H) = 10.0, log i? = 18.0). The contour plot shows that in the low density 
limit {logn{H) ^ 4.5 cm~^) the ratio is a good temperature indicator in that moving 
perpendicular to the ionization parameter the ratio changes are roughly independent of 
density. At higher densities the ratio becomes dependent upon both temperature and 
density. 

The gas density integrated line ratio as a function of distance shows that for large 
distances from the ionizing source the ratio is independent of the weighting of the density 
axis; i.e., the ratio is in the low density limit. Using the analytical expression from 
Osterbrock (1989) the ratio can be converted into temperature in the low density limit. 
For the range 20.0 cm < logi? < 22.0 cm the inferred temperature has the range 14,000 K 
> Te > 8400 K, falling off very gradually as approximately R"^^^. At smaller radii the 
A5007/A4363 ratio further decreases, reflecting, mainly, the drop in A5007 due to coUisional 
quenching and the increase of A4363 which continues to emit efficiently up to its critical 
density, rather than reflecting an increase in temperature. 
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The [Nil] A5755/A6584 line ratio is also a useful temperature indicator. Wilson, 
Binette k, Storclii-Bergmann (1997) use this ratio in concert with the [Olll] A5007/A4363 
ratio as evidence for the presence of matter-bounded clouds in the NLR. Defining R[nii] 
as the [N ll] line ratio and R[oiii] as the [O ill] line ratio, our model integrations predict 
that R[oiii] versus R[nii] depends little upon the density weighting factor g{n). We find a 
linear-log relationship of the form 

logi?[0777] = 0.1593 log i?[^77] + 0.1932. (3) 

The fiducial dust-frcc integration from Figure 3 (§3.1) predicts R[oiii] = 0.0182 and 
RlNii] = 0.0215, and the fiducial dusty simulation predicts R[oiii] = 0.0206 and 
R[Nii] = 0.0229. These correspond to temperature differences {Tpiii] — T[nii]) of 3000 K 
and 3500 K, respectively. These predictions lie in the range of observed values (Wilson, 
Binette & Storchi-Bergmann) , and do not require the presence of a significant population 
of matter-bounded clouds, although we do not exclude this possibility. 

5.3. Ionization 

The ratio [Olll] A5007/(Hq; -|- [Nil] A6584) has been used as an ionization indicator 
and this ratio is shown in the top left panel of Figure 5b. The outer contour has the 
linear value of 1.0 and the dashed fines are -1-0.2 dex steps. We have chosen to include 
both Ha and [Nil], since many filters include both lines. The contour plot of this ratio is 
relatively flat over a wide range in density, distance and ionization parameter. In fact, for 
a constant density of n{H) = cm^^, the [OlIl]/(Ha + [Nil]) ratio would be constant 
from ~ 10 - 50 1%^^ pc, even though the ionization parameter decreases by ~ 10 over 
the same range. A ratio that behaves in similar manner is [O ill] A5007/H/5, and is shown 
in the bottom left panel of Figure 5b. The outer contour has the linear value of 1.0 and 
the dashed lines are -1-0.2 dex steps. The [Olll]/H/3 ratio is just as fiat over a wider range 
in ionization parameter as is [O ill] /(Ha -|- [Nil]). Thus, while these fine ratios are useful 
in separating Hll region galaxies from active galaxies (c.f. BPT), they are not reliable 
ionization parameter indicators. 

In the top right panel of Figure 5b the [Olll] A5007/(Hq; + [Nil] A6584) as a function 
of distance is shown. The ratio has the feature of being constant from ~ 3 - 300 L^^^^ pc 
for the integration assuming a steep distribution {g{n) oc n~^^^■, dashed line) of clouds along 
the density axis. The solid line {g{n) oc n~^ °) is flat over a much smaller distance range, 
because in this integration emission from higher density gas is being included. The sharp 

1 /2 

falloff of the ratio at large distances {R > 300 I/43.5 pc) occurs because the [O ill] A5007 
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emissivity is dropping very quickly with distance at the lowest gas density (10^ cm~^) 
considered here, since the gas is underionizcd and can not produce the line efficiently, while 
the [Nil] A6584 line is simultaneously increasing in strength (see Figure 1). The bottom 
right panel of Figure 5b shows the integrated [Olll]/H/3 ratio as a function of distance; 
this ratio behaves similarly to the [OlIl]/(Ha; + [Nil]) ratio. Both of these line ratios are 
observed to be relatively constant over factors of 3 - 4 in distance in several Seyfert 2 
galaxies. (Robinson et al. 1994; Capetti et al. 1996). 

Figure 6a shows the [Olll] A5007/[Oll] A3727 line ratio as a contour plot. In the 
figure, the outer contour has a linear value of 1.0, and the dashed lines arc +0.2 dcx steps. 
At low densities, the contours are closely spaced and parallel with ionization parameter, 
indicating that this hne ratio is a good ionization indicator. At densities ~ 10^-^ cm~^, [On] 
becomes coUisionally quenched, while [Olll] does not, thus the ratio becomes dependent 
upon density as well. 

In Figure 6b we show the predicted [O ill] /[Oil] line ratio as a function of radius for 
4 models compared with the observations of Robinson et al. (1994) for NGC 4151. The 
data have been placed on the distance axis by assuming that 1".0 ^ 100 pc, as used by 
Robinson et al. The model integrations have been rescaled in distance by factors of 4, 
corresponding to a 16-fold increase in the ionizing luminosity incident upon the clouds. 
The observed strengths of these and a number of other emission lines in NGC 4151 suggest 
that the ionizing luminosity incident upon the narrow-line clouds may be 10 - 20 times 
the observed one (Penston et al. 1990; Robinson et al.), though the magnitude of the 
photon deficit is very uncertain (see also Schulz & Komossa 1993). Since the ionizing 
luminosity employed here happens to match that of NGC 4151 emitted toward Earth 
(Robinson et al.), our rescaling is consistent with the continuum beaming inferred from the 
observations. The predicted [O ill] /[Oil] ratio shown in Figure 6b did not include emission 
from matter-bounded clouds; including such clouds would tend to enhance the [O ill] /[Oil] 
ratio where it is small in Figure 6b. 

Robinson et al. (1994) include several other line ratios as functions of projected distance 
for NGC 4151. These include [Ne v]/[Nelll] and [O ill] /[Nil] which behave similarly to 
[Olll]/[Oll] shown in Figure 6b. Also observed are [Nii]/Hq;, [Ol]/Ha, [Oii]/Hq;, and 
[Sii]/Hq;. Our simulations qualitatively match the general trends in these line ratios. 
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6. Is the NLR Dusty? 

The presence of dust in the NLR is a vital clue to the origin of the gas. For instance, if 
shocks are prevalent in the emitting region, then grains will not survive (Donahue & Voit 
1993). Both the dust-free and dusty simulations reproduce the average narrow emission line 
ratios shown in Figure 3. However, other lines do indicate substantial depletions relative to 
solar. The absence of an observed [Call] A7291 line has led to speculation that the NLR is a 
dusty place (Ferland 1993; and Kingdon, Ferland & Feibelman, 1995). Recent observations 
of the [Call] line (Villar-Martm & Binette 1997) show that the line is much weaker than 
the solar integrations predict (Table 1), supporting Villar-Martin & Binette conclusion that 
Ca is sharply depleted. The simulations presented here indicate that Ca is depleted relative 
to solar by factors of 3 - 160 based on the range of observed line strengths when compared 
to the dust-free simulations in Table 1 . The results of the simulations for the infrared [Fe ll] 
lines, as discussed in §3.1.3, indicate that Fe is also underabundant relative to solar by 
factors of 10 - 20 in NGC 4151. Comparisons of the compilations of the [Fell] 1.257/xm/P/? 
line ratio by Simpson et al. (1996) and Veilleux et al. (1997) with the simulations provide 
further evidence that the NLR is a dusty environment. A similar conclusion was reached by 
Simpson et al. 

Two line ratios, La/HP and Mgll/H/3, observed in NGC 3393 and discussed in §3.1.1, 
also suggest that dust exists in the emitting gas. However, the uncertainty in the correction 
for the unintended positional error in the observations, makes the conclusion from Lq;/H/3 
uncertain. There is no positional error for the Mgll/H/3 ratio, and this ratio makes a 
stronger statement supporting the presence of dust in NGC 3393. 

Evidence that a portion of the emitting gas is dust free is apparent in the observations 
of strong coronal lines in many AGN. Oliva et al. (1994) have observed a very strong 
[Cavill] 2.32 ^m infrared line and strong optical [Fevil] A6087, [Fex] A6374, [Fexi] A7892, 
and [Fexiv] A5303 lines, indicating that the gas is dust-free. Simulations of coronal line 
emitting gas by Korista & Ferland (1989), Oliva et al. (1994) Ferguson, Korista, & Ferland 
(1997), and the results of those shown in Table 1 indicate that these elements are not 
depleted. 

Taking the observations at face value, we are presented with a conundrum: the low 
ionization UV Mgll, optical [Call], and infrared [Fell] lines clearly indicate that dust 
is present, but the strong coronal lines suggest a dust-free environment. One possible 
reconciliation is to postulate that there is dust in the NLR, but only in the neutral or 
partially neutral portions of the clouds. In the simulations presented in this paper we have 
assumed that the entire cloud is either dust-free or dusty. It is not unreasonable to imagine 
that the highly ionized portion of NLR clouds to be dust free, perhaps because the grains 
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have been exposed to a strong UV continuum for quite some time, and that in more neutral 
regions grains are free to exist. Another possibihty is that high ionization clouds have 
been processed through shocks, but those which emit more neutral ions have not. Both of 
these possiblities argue for a very complex NLR environment, where the emitting gas has a 
variety of origins. 

7. Summciry 

We have shown that the AGN narrow emission line spectra, from [Ol] to [Siix], can 
largely be reproduced by assuming that the emitting region consists of clouds with a 
wide range of gas densities and distances from the ionizing source. Here we chose simple 
power-law cloud distribution functions and the appropriate power law indices {(3 and 7) 
whose resulting spectrum matched the observed [Olll]/[Oll] and [Olll]/H/9 ratios. While 
smooth power law distributions in cloud properties may be too simple, its success in 
predicting the observed spectrum points to the strength of the "locally optimally emitting 
cloud" scenario. The LOG model integrations result in the following: 

(1) The classical BPT line ratios for AGN narrow emission lines are reproduced (Figure 3). 

(2) The predicted H/5 equivalent width, assuming full source coverage and an isotropic 
continuum, is more than sufficient to account for the observed emission in most Seyfert 
galaxies (Table 1). An anisotropic continuum will strengthen the predicted line equivalent 
widths and allow for smaller covering fractions. The predicted [Qui] /[Oil] ratio as a 
function of radius implies that the narrow line clouds see a continuum that is an order of 
magnitude more luminous than we observe in NGG 4151. 

(3) The LOG integration simultaneously reproduces the low [Sll] A6716/A6731 density 
and the high [Olll] A5007/A4363 temperature (cf. Figure 3b), resolving a long standing 
problem. 

(4) The observed line width - critical density correlation for the optical forbidden lines 
can be reproduced. A natural consequence of an ensemble of clouds in the gas density - 
distance plane is that lines with higher critical densities are formed at smaller distances (cf. 
Figure 1). We have shown that if we assume that this ensemble of emitting clouds is placed 
in a gravitational potential, then the line width - critical density correlation is very similar 
to the observed will naturally result. 

(5) The presence of a large density range at each radius can change the values of the 
classical density and ionization-level indicators to at least more closely mimic the observed 
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radial gradients in these intensity ratios. 

(6) Low ionization lines show clear evidence of depletion (Figure 3 and Table 1), but high 
ionization regions do not. It is possible that dust grains are selectively destroyed in high 
ionization regions by shocks or sublimation and that the neutral regions are undisturbed. 
This argues for a variety of NLR cloud origins. 

We thank Andrew Cooke for agreeing to our use of the NGC 3393 data. We acknowledge 
suport from NASA and NSF through NAG-3223 and AST93- 19034, and STSCI for GO- 
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its director John Connolly for access to the Exemplar parallel supercomputer. 
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Figure Captions 

Fig. 1. — Contours of constant logarithmic line equivalent widths as a function of logR and 
\ogn(H) for the 23 ions indicated, referenced to the incident continuum at 4860 A. The 
bold lines represent 1 dex increments and the dotted lines are 0.2 dex steps. The triangle is 
the peak of the equivalent width distribution and the contours decrease downward to the 
outer value of 1 A. The reader will sometimes find it convenient to view the contour plots 
along the ridge at large inclination angle to the sheet of paper. The upper right hand plot 
on the first panel of the figure is the log Tg of the front face of the cloud. The temperature 
decreases from 10^ K in the lower left hand corner of the plot to 10^ K in the upper right 
hand corner. Again, the bold lines represent 1 dex increments and the dotted lines are 0.2 
dex steps. 

Fig. 2. — Same as Fig. 1, except for the dusty abundance set described §2.2.2. 

Fig. 3. — Line diagnostic diagrams after the manner of Baldwin, Phillips & Tcrlcvich (1981). 
Difi^erent line types represent different abundance types, dashed are solar or dust-frcc, and 
solid lines are the dusty abundance set described in the text. Different lines represent 
different families of the density weighting index, /? with different covering fraction indices, 
7 along each line with the arrows indicating the direction of the more negative index. The 
star represents the same model integration of dusty abundances with indices /(r) oc r~^^^ 
and g{n) oc n~^'^. The square represents the same model integration with solar abundances 
with indices /(r) cx r^^'^^ and g(n) oc n~^'^. The large "S" and "L" symbols refer to the 
observations of the mean Seyfert type 2 objects and LINER type objects, respectively, as 
discussed in §3.1. 

Fig. 4. — Line width versus critical density of the square point shown in Fig. 3. The solid 
circles are the ftill width at 10% intensity with a fitted slope of 0.193. The open squares are 
FWHM with a fitted slope of 0.118. The observed range in the slope is from 0.095 to 0.200. 

Fig. 5. — Line ratio diagnostic diagrams of gas density, temperature, and ionization shown 
in two different forms. The left hand column of plots in the figure are logarithmic contour 
diagrams of the ratios indicated taken from the simulations shown in Figure 1. The right 
hand column of plots in Figure 5 are graphs of the indicated line ratio versus distance from 
the central source for two different integrations as discussed in §3, but integrated in the 
density dimension only. The solid lines are integrations assuming that g{n) oc n^^'^, and 
the dashed lines assume g{n) oc n''^^^. Each line ratio is shown on a different scale and we 
will discuss each in turn in the text in §5. 
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Fig. 6 — The predicted [Olll] A5007/[Oll] A3727 ratio shown in two forms: (a) is a 
logarithmic contour diagram of the ratio taken from the simulations shown in Figure 1. 
The outer contour has a linear value of 1, and the solid lines are 1.0 dex increments with 
the dashed lines being 0.2 dex steps, (b) The line ratio as function of distance from the 
ionizing source. The lower solid line is an integration in density assuming g{n) oc n~^'°, and 
the upper solid line assumes g{n) oc n~^'^. The dashed line is the fiducial integration with 
g{n) oc n"^-^ for the dust- free simulation. The dotted line is the fiducial integration with 
g{ri) oc n"^'^ for the dusty simulation. All integrations have been rescaled by factors of 4 in 
distance corresponding to 16 in luminosity (see §5.3). The filled circles are the observations 
of NGC 4151 by Robinson et al. (1994), assuming 1".0 ^ 100 pc. 
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TABLE 1 



Observed versus integrated spectra 



Line 


Mean Sy 2« 


Reference'' 


NGC 3393^ 


Circinus'^ 


Solar^ 


dusty-'' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Wa(H/3)^ 


5-30 


N 


— 


— 


152 


61 




15 - 100 


BFP 


— 


— 


— 


— 


La A1216 


55 ± 20 


FO 


20.9 


— 


28.7 


9.62 


Nv A1240 


— 


— 


1.42 


— 


0.43 


0.60 


Si IV A1397 


— 


— 


<0.58 


— 


0.09 


0.03 


Oiv] A1402 


— 


— 


<0.84 


— 


0.22 


0.45 


C IVA1549 


12 ± 8 


FO 


4.75 


— 


3.70 


2.47 


Hen A1640 


— 


— 


1.75 


— 


2.18 


2.82 


Cm A1909 


5.5 ± 3.7 


FO 


1.60 


— 


1.25 


3.63 


Cii A2326 


- 


- 


0.25 


- 


0.39 


0.97 


[Neiv] A2424 






0.32 




0.19 


0.31 


[0 II] A2471 






0.11 




0.16 


0.17 


Mgii A2798 


1.8 ± 1.5 


FO 


0.43 




1.42 


0.23 


Hen A3205 






0.13 




0.12 


0.15 


[Nev] A3426 


1.2 ± 0.2 


FO 


0.36 




0.62 


0.67 


[On] A3727 


3.2 ± 1.7 


FO 


2.27 




2.32 


3.20 


[Neiii] A3869 


1.4 ± 0.4 


FO 


0.81 




1.13 


1.00 


[Sii] A4074 


0.3 ± 0.2 


FO 


0.12 




0.26 


0.28 


m A4102 






0.18 




0.27 


0.27 


H7 A4340 






0.48 




0.48 


0.47 
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TABLE 1 Continued 
Observed versus integrated spectra 



Line 
(1) 


Mean Sy 2« 
(2) 


Reference'' 
(3) 


NGC 3393^ 
(4) 


Circinus*^ 
(5) 


Solar^ 
(6) 


dusty-'' 
(7) 


[Oiii] A4363 


0.2 ± 0.1 


FO 


0.09 




0.21 


0.22 


He II A4686 


0.3 ± 0.1 


FO 


0.25 


0.30 


0.30 


0.37 


Rp A4861 


1.00 


FO 


1.00 


1.00 


1.00 


1.00 


[Oiii] A5007 


10.8 ± 3.0 


FO 


10.2 


11.8 


11.7 


10.8 


[Ni] A5199 


0.15 ± 0.09 


FO 


0.15 


0.35 


0.07 


0.28 


[Fexiv] A5303 








<0.10 


0.033 


0.004 


He I A5876 


0.13 ± 0.06 


FO 


<0.11 




0.13 


0.11 


[Fevii] A6087 


0.10 ± 0.05 


FO 




0.055 


0.049 


0.009 


[Oi] A6300 


0.57 ± 0.20 


FO 


0.35 


0.36 


0.56 


0.84 


[Fex] A6374 


0.04 ± 0.04 


FO 




0.010 


0.015 


0.002 


Ha A6563 


3.1 ± 0.1 


FO 


3.00 


2.73 


2.96 


2.91 


[Nil] A6584 


2.9 ± 1.0 


FO 


3.64 


3.19 


1.18 


1.69 


[Sii] A6720 


1.5 ± 0.5 


FO 




1.64 


1.20 


1.87 


[Arm] A7135 


0.24 ± 0.07 


FO 




0.16 


0.22 


0.27 


[Can] A7291 


0.04 ± 0.04 


VB 






0.16 


0.002 


[On] A7325 


0.10 - 0.46'^ 


DS 




0.073 


0.20 


0.21 


[Fexi] A7892 








0.066 


0.068 


0.008 


[Siii] AA9069,9532 


0.5 - 4.5'^ 


DS 




0.94 


1.40 


1.61 


[Sviii] A9915 








0.044 


0.024 


0.022 


[Six] 1.252 








0.042 


0.11 


0.086 


[Fell] 1.257 /im 


0.12 ± 0.01 


T 




0.09 


1.16 


0.23 


[Fell] 1.279 /ill! 


0.03 ± 0.01 


T 






0.05 


0.01 


Pp 1.2818 urn 


0.16 ± o.or 


T 




0.15 


0.16 


0.15 
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TABLE 1 Continued 
Observed versus integrated spectra 



Line 


Mean Sy 2° 


Reference'' 


NGC 3393^ 


Circinus*^ 


Solar^ 


dusty-^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


[Fell] 1.534 //m 


<0.02 ± 0.01 


T 






0.10 


0.01 


[Fell] 1.644 jjm 


0.09 ± 0.01 


T 




0.07 


1.12 


0.23 


[Fell] 1.677 /iin 


0.02 ± 0.01 


T 






0.07 


0.01 


[Sivi] 1.96 fim 








0.088 


0.088 


0.015 


Bvy 2.166 /xm 








0.027^ 


0.029 


0.026 


[Cavill] 2.32 /xm 








0.065 


0.15 


0.002 


[Sivil] 2.48 /xm 








0.15 


0.11 


0.020 


[Siix] 3.94 urn 








0.19 


0.17 


0.027 


Bra 4.052 pim 








0.079^ 


0.080 


0.073 


[Fell] 25.98 pim 










0.64 


0.13 



"The range of observed values taken from the reference in column 3. 

''reference of the observed values: (N) is Netzer 1990, (BFP) is Binette, Fosbury, & Parker (1993), (FO) 

is Fcrland & Osterbrock (1986), (VB) is from Villar-Martm & Bincttc (1997), (DS) denotes a range of 
values from Dopita & Sutherland (1995), and (T) is from an observation of NGC 4151 by Thompson 
(1995). 

'^Observations of NGC 3393 from Cooke et al. (1997). 
''Observations of the Circinus galaxy from Oliva ct al. (1994). 

'^Predictions from the integrated spectrum assuming solar abundances indicated by the square in Fig. 3. 
■''Predictions from the integrated spectrum assuming Hll region abundances indicated by the star in 
Fig. 3. 

^Equivalent width of H/3 in A. 

''These observations are given as a range. 

'Assuming Case B conditions. 
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